RESONANCE SPECTRUM FOR ONE-DIMENSIONAL LAYERED MEDIA 

ALEXEI lANTCHENKO 

>—^ 1 Abstract. We consider the "weighted" operator Pk = —dxa{x)dx on the hnc with a step- 

CN ' like coefficient which appears when propagation of waves thorough a finite slab of a periodic 

►^,^^' medium is studied. The medium is transparent at certain resonant frequencies which are 

frt , related to the complex resonance spectrum of P^ ■ 

If the coefficient is periodic on a finite interval (locally periodic) with k identical cells then 
the resonance spectrum of Pf; has band structure. In the present paper we study a transition 
to semi-infinite medium by taking the limit fe — >■ oo. The bands of resonances in the complex 
lower half plane are localized below the band spectrum of the corresponding periodic problem 
(fe = oo) with fc — 1 or fc resonances in each band. We prove that as fc — >■ oo the resonance 
spectrum converges to the real axis. 



\o 



'^ • 1. Introduction 

In the present paper we consider operator Pk — —dxak{x)dx on the line with step-hke coeffi- 
cient Ofc which is periodic on a finite interval defined as follows: 

>(1) ak{x) = a{x), for x G [0, fc]; ak{x) = --^, for x ^ [0, fc], 

^^ , where a{x) is 1— periodic function equal to 

^: (2) ao{x)^i % forxejO,..) 

^^ \ for a; e [0, 1). Here &i^2 > and < a;2 < 1- Equation 

Pki^ = -dxak(x)dxi^{x) = A^-0 

appears when the propagation of waves through a finite slab of a periodic medium is studied. 
Jv>( \ Such systems are also called finite or locally periodic media (for revue see 6 ). 

}_( ' When fc is large then the properties of medium is close to an infinite periodic problem in a 

5t ! sense that we are going to discuss in the present paper. 

We denote P = —dxa{x)dx the pure periodic operator, where a{x) is 1— periodic function 
equal to flo for x E [0, 1) as in ©. Then the Floquet theory shows the existence of a pair of the 
quasi-periodic solutions ip± of the equation —dxa{x)dxip± — A^V'±j 

V'^(A,a; + l) = e±^V=^(A,a;), 

such that ip^ e i^(K±) for ImA > 0. Here 9 = 9{X) is the Bloch phase. We denote 

(3) F(A) = ^ cos{A(a;2fo2 + (1 - .t2)&i)} - ^ cos{A(x252 - (1 - 2:2)61)} 
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the Lyaponov function for P (see Section [^TT|) . Here 

bl + bl 
P- 



2&162 
The spectrum of the operator P has band structure with allowed zones defined as follows; 

(4) Xea{P) ^ \F{\)f < 1, AeM 

(see [4] and Section 14]). The band edges are given by solutions of F{X) = ±1. 

The relation between the Bloch phase 6 and the spectral parameter A is called dispersion 
relation: 

cos6'(A) ^F{X). 

Since the coefficient a{x) is constant equal to l/b^ outside a finite region, we are here concerned 
with a scattering problem. 

We shall denote the reflection and transmission coefficients for the operator P^ by r^ and tk, 
respectively: 

PfcV = AV, V- = e'^'"^ + rfce-*^''l^ x < 0; V = ifce'^^'i", x > k. 

Following the ideas in [7] we consider a transition to semi-infinite periodic materials by taking 
the limit fc — ;> cxo of the reflection coefficient r^ for P^ . 

The limiting operator 

Poo^ = ~dxaoo{x)dx^{x) 
corresponds to the case of such a long slab that it can be considered as half infinite. 

In the case of the operator P.^, the solution %p of the scattering problem is defined as the 
solution of the equation PocV' = -^^V-"! such that 

(5) i) = e*^''!^ + re-'^^^"", x < 0; ?A = cV'+(A, x), x > 0, 

with some r = r(A), c = c(A). 

As in [7] we have that the reflection coefficients rfc(A) and r(A) are analytic in the upper half 
plane C+ = {A : ImA > 0} and continuous in C+, and rfc(A) — > r(A) when k ^> 00 and A e C+. 
When A is real, r/c(A) converges to r(A) in the weak sense (see Theorem[2l Section [Z!2|) . 

Numerical calculation shows that in each allowed zone of (j{P) there are in general k — 1 
frequencies Xj where the transmission probability is one: \tk{Xj)\^ = 1, j' = I, . . . , fc — 1, and the 
medium is perfectly transparent: |?'fc(Aj)p — 0. There exist an additional frequency Aq when the 
medium consisting of only one unit cell is transparent and then |tn(Ao)P = 1 for all n = 1, 2, . . . , A;. 
The pics in the transmission probability are related to the complex resonances close to the real 
axis. 

We make the following definition. 

The operator P^ defined from {u e H^{^), ad,^u e H'^iR)} to L^{R) is self-adjoint. For 
ImA > 0, we call Rk{X)v = {—dxa{x)dx — X^)^^v the resolvent of Pk- For any k — I, 2, ... , the 
operator-valued function 

Rk{X) : Ll,^^{R) ^ LlM) 
can be continued to the lower complex half-plane C_ as a meromorphic function of A G C and it 
has no poles for ImA > — efe, A 7^ 0, with efc > positive constant dependent on k (see Section 

n. 

The poles of the i?^ (A) in C_ are called resonances or scattering poles. We denote the set of 
resonances Res [Pk)- 

Using the explicit construction of the resolvent in [4] the poles are calculated numerically. 
Some examples are presented in Section [^ figures ([T]), ([2]) and ([3]). We summarize the properties 
of Res {Pk) in the following Theorem. 
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Theorem 1. We consider the finite periodic Pk and periodic P operators generated by the same 
unit cell given in 0^. Let Res (Pk) C C_ denote the resonance spectrum for the finite periodic 
system, with k — 2,3... identical cells and <j{P) C M denote the band spectrum for P given by 

The resonance spectrum for the finitely periodic system has band structure related to the bands 
of the real spectrum for the pure periodic problem as follows: 

1) The resonance spectrum of Pk has band structure. Resonances are localized below the bands 
of the real spectrum of P : 

AeRes(Pfe) ^ Re (X) satisfies 1^ ^ Re (A) G cr(Poc). 

Each resonance band of Pk consists ofk—1 resonances Ai, . . . , A^ — 1 and eventually an additional 
resonance with real part ReA = Aq = Trm/x2b2, to € Z, such that |ii(Ao)P = 1, the one-cell 
medium is "perfectly transparent" at frequency Aq. 

2) If the condition 

(6) b2X2 ^ bi{l - x^) ^ -^ = — ^. 

X2 1-2^2 

is satisfied then Aq = TTm/x2b2, m (z 1j, is the degenerate band edge (two bands has common edge 
at Xq). The resonance spectrum Res (Pk), k — 1,2,3, ... , is periodic with the period T = t^^- 

3) As /c — )■ oo then the resonance spectrum of Pk approaches the real axis. 

In the present paper we motivate these numerical resuhs. 

The band structure of the resonance spectrum for a finitely periodic system and its relation 
to the band spectrum of the correspondent periodic problem is well-known in physical literature 
(see n). 

We say that Res {Pk) is periodic if there exists T > 0, period, such that 

Res (Pfe) n {[q + Tn,p + Tn] - iM.) = Res {Pk) n {[q + Tm,p + Tm] - iM) 

for any q < p and n,m € Z. This property follows directly from the equations defining the 
resonances in Section (|4.2I) if condition ^ is satisfied. 

A special property of the operator P = —dxa{x)dx with step-like periodic coefficient a{x) is 
that the coefficients (r ± 1) in the dispersion relation 

(7) 2 cos6i(A) = (p + 1) cos{A(j:252 + (1 - 2:2)61)} - (p - 1) cos{A(2;262 - (1 - X2)bi)} 

are independent of the spectral parameter A. Formula ([7]) implies that the band spectrum is 
periodic if the profile of a verifies ^. 

The third part of the Theorem is proved in Section [6] 

The convergence of the resonances for a finitely periodic system with k cells to the bands of 
real spectrum for the periodic problem as /c — >■ cxd was discussed by F. Barra and P. Gaspard in 
[3] in the case of Schrodinger equation. 

In our proof we use representations for the reflection and transmission coefficients rk , tk for 
a finite slab of periodic medium as in the recent paper of Molchanov and Vainberg [7]. The 
authors considered transition of truncated medium described by the 1— D Schrodinger operator 
to semi-infinite periodic materials. By relating the reflection coefficients to the resolvent of Pk 
we show explicitly that the resonances correspond to the poles of the analytic continuation of 
rk{X) to C_. Then we consider the limit of the poles of rk{X), as fc — >■ 00. 

Note that for A S R the reflection coefficient rk+i for fc + 1 cells medium is related to r^ for 
k cells medium via rk+i = f\{rk), where f\ is a linear- fractional automorphism of the unit disk. 
By considering the fixed point of /> we get a new proof of the convergence of rk{X) to r(A) when 
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A belongs to the spectral gapes and non-degenerate band edges for the operator P (see Section 

m- 

The structure of the paper is the following: 
In Section [2] we recall some well-known facts concerning spectral problem for weighted Sturm- 
Liouville operators (see [2]) and consider scattering by a finite slab of a periodic medium. We 
follow [7] with minor changes due to the special form of operator Pk . We recall exact formulas 
for the the reflection and transmission coefficients using the iteration of the monodromy matrix. 
We recall also the result of [7] on a transition to semi- infinite periodic material (limit k — > oo). 
In Section [3] we give explicit expression for the monodromy matrix of Pk ■ In Section [5] we recall 
the iterative procedure used in j4: for construction of the resolvent Rk and define resonances. 
In Section [5] the reflection coefiicient r^ is expressed using the iteration formulas of [? and we 
show that the poles of i?fc(A) and the poles of rfe(A), A G C_, coincide. In Section [HI we prove 
the convergence of Res(Pfe) to the real axis. In Section [7] we discuss the convergence of rfe(A), 
fc -> oo, for A G M by considering the limit of a sequence of linear- fractional automorphisms on 
the unit disk. In Appendix lAl we present numerical examples. 

Acknowledgements. The author would like to thank Maciej Zworski for suggesting to look 
at the problem considered in the present paper and for helpful discussions. 

2. General methods for truncated periodic operators 

In this section we following 7 consider the scattering theory for operator Pk combining the 
Floquet-Bloch theory and scattering theory for 1— D weighted operators. 

2.1. The monodromy matrix and Bloch quasi-momentum. We recall first some well- 
known facts concerning the spectral problem of Sturm-Liouville operators on the line (see [2]). 
We consider equation 

(8) PV = -dMx)dri^ix) = A^V' 

on {tp G iJ[Q^(M), adxip G iJ[Q^(M)} with a strictly positive a{x) as in the Introduction, formula 
((H) or periodic as qq in ^. 

Let -01,2 be solutions of ([8]) with initial data 

(9) V'i(A,0) = 1, (aa,0i)(A,O) = 0; M\0) = 0, (a9,i^2)(A, 0) = 1. 
We define the transfer matrix (propagator) Mx{Q,x) for operator P 

(10) M,{Q,x) = (^ (aa.^.)(A,.) („5^^^)(;,^^) j • 

From ^ it follows that M\{{), x) is the identity matrix. For any solution -0 of (|5]) matrix Ma(0, x) 
maps the Cauchy data of at x = into the Cauchy data of ■0 at point x : 






A / V A 

As the generalized Wronskian associated with ^ji , 02 

I^[V'l,02] = i>iadx^2 - '02a9a;0i 

is constant, we have 

detMA(0,x) = VF[0i,02](l) = M^[0i,02](O) = 1. 

Equation ((S)) with ImA > has exactly one solution t/A in L^(R+) normalized by the condition 
0"'"(A, 0) ~ 1, and it has exactly one solution 0^ G i^(M_) normalized by the same condition. 
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Here R± are semiaxes x ^ 0. Any solution of (|S]) can be represented as linear combinations of 
tpi and "02 and from the normalization of ip it follows, that there exist functions m — 77i^(A) 
such that 

■0± =-01 +m=^(A)V'2, ImA > 0. 
Functions m^ — m^(A) are called Weyl's functions and we have 

V'i+m+(A)02 e i2(M+), 01 +TO"(A)V'2 e i^(K_), ImA>0. 

Let a{x) be periodic: a{x + I) — a{x). Consider propagator through one period (monodromy 
matrix) : 

Af,-M,(0,/)-(^ ^ ^ j (A) - (^ (oa^V^KM (a9,V2)(A,0 j " 

Denote F{X) = i tril/A = {a + S)/2 = (■0i(A, 1) + aa02(A, 1) the Lyapunov function. Both Mx 
and F{X) are entire function of A and det Mx = 1. The eigenvalues A'^(A) of Mx are the roots of 
the characteristic equation 

(11) M^ - 2^F(A) + 1 = 0. 

If ImA > then one can select roots fi^{X) of pip in such a way that /i^(A) — e^*'^^'^', where 
9{X) is analytic and 

(12) Im6'(A) > when ImA > 0, 
i.e., 

(13) Im+(A)| < 1, |/^"(A)| > 1, ImA > 0. 

The roots /i (A) for real A > are defined by continuity in the upper half plane: 

Ai^(A) =//±(A + iO), Ae [0,oo). 
Since the trace of Mx is equal to the sum of the eigenvalues e^'^^^^^\ 

(14) cosW^F{X) = i(0i+aV'2)(A,/) = ^(a + (5). 

The spectrum of P belongs to the positive part of the energy axis E — X^ and has band 
structure. 

For real A > 0, the inequality |-F'(A)| < 1 defines the spectral bands (zones) 

&n = [A2n-1, A2„], 71=1,2,..., 



on the frequency axis A = vE. The bands are defined by the condition |F(A)| < 1 and F{X) = ±1 
at any band edge A = Aj . 

The function 9{X) is real valued when A belongs to a band. The roots /i^(A) are complex 
adjoint there, and |/i^| = 1. 

The spectrum of P (on L'^{M)) on the frequency axis is U^i ^n- 

The complimentary open set, given by |-F(A)| > 1, corresponds to spectral gaps, IJ^i 9n- On 
gaps, the function id{X) is real valued, the roots fi^{X) are real and ([T3l) holds. 

A point Xj which belongs to the boundary of a band and the boundary of a gap is called a 
non-degenerate band edge. If it belongs to the boundary of two different bands, it is called a 
degenerate band edge. 

As in [7] we get that if A = Aq is a non-degenerate band edge, then F'{X) 7^ 0. If A = Ao is a 
degenerate edge, then F'{X) = 0, F"{X) 7^ 0. Both eigenvalues of the monodromy matrix AIx at 
any band edge are equal to 1 or both are equal to —1. 



h^{\) ^ , _±, 
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We normalize the eigenvectors h^{X) of Mx by choosing the first coordinate of h^{X) to be 
equal to one: 

1 

m±(A) 

The second coordinates of the vectors h^{X) coincide with the Weyl's functions defined above. 
In fact, if ip^ are solutions of the equation P^|J ~ X^ip with the initial Cauchy data given by the 
eigenvector /i^, then 

(15) V^(A,x + 0=e±''''(^V*(A,x) 

and (fT5|) implies that -0^ G L?(J^±) when ImA > 0. From here it follows that tj}^ coincide with 

Weyl's solution introduced for general Hamiltonians P, and that the second coordinates of the 

vectors /i^(A) are Weyl's functions. 

Since 

a - e±^'^(^) p \ / 1 

the following two representations are valid for Weyl's functions: 

(16) m±(A) = ^ > - ^^ ' 



0, 



2.2. Reflection coefficient for the truncated periodic operator. We consider operator Pk 
with the truncated periodic coefficient au '■ 

Pkip = ~dxak{x)dxtpix), ak{x) = a{x), for x G [0, fc/]; ak{x) = -pr, for x ^ [0, kl], 

H 

which appears when the propagation of waves through a finite slab of a periodic medium is 
studied. We shall also consider the limiting case k — oo : 

-Poo'0 = -dxaooix)dx^p{x), 

which corresponds to the case of such a long slab that it can be considered as half infinite. 

We shall denote the reflection and transmission coefficients for the operator Pk (with compactly 
supported coefficient a^) by r^ and ife, respectively: 

PfcV = AV, V' = e'^''^'-' + rfee-*^''^^, x < 0; V = ifce*^''^^ x > kl. 
In the case of the operator Poo, the solution ip of the scattering problem is defined as the 
solution of the equation Poo'0 = A^f/'j such that 

(17) ^ ^ e'^'"'' + re-'^'"'', x < 0; ?A = cV'+(A, x), x > 0, 

with some r = r(A), c — c(A). We have the following version of Theorem 3 of S. Molchanov, B. 
Vainberg in [7]: 

Theorem 2. 1) The transfer matrix over k periods M^ = T\{0, Ik) has the form 

.1„^ ,.k ( «fc Pk \ sinfc/g(A) sin(fc - l)ie{\) 

where 9 — 0{X) is the Bloch function. The elements of M^ satisfy the relations 

smkW{X), ,, ^ sinfc/6'(A)^ 

ak-6k^ . , - (a~d , Pn = . ,„.,. P, 
sm/6'(A) sm/t^(A) 

, . sinfc/0(A) , , 

19 7fe = — Wrr7, «fe + 4 = 2cosfc;0 A . 

sm/t^(A) 
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2) The reflection coefficients have the forms 

(a-<5)+*(6i7+^) 



(20) rfc(A) = - 



_ -^+bij-i{a-S) 



(21) r{X) - „ 

2sin/0(A)-(6i7-^: 

3) The transmission probability have the form 

(22) M\)\'^ 






1 



4) r/ie reflection coefficients rfc(A) and r(A) are analytic in the upper half plane C+ = {A : Im A > 
0} and continuous in C+. For any A € C+ \ U5^]^6„ we have rfc(A) — >■ r(A). When A G U^]^6„, 
rfc(A) converges to r(A) m i/ie iweafc sense: 



rfc(A)(p(A)(iA — > / r(A)iy9(A)(iA as fe — > oo. 

; -/ — oo 

/or any iesi function Lp ^ D. 

Proof: We reproduce here the proof of 7 for the sake of completeness with only minor 
changes due to the "weight" in the definitions of Pk and Poo ■ Formula (|18p follows by induction 
from relation pT|) : 

Ml - 2 cos W{\)Mx +1 = 0. 

The first three relations of ([T9|) are immediate consequences of ([T8| . In order to get the fourth 
one we note that the eigenvalues of M^ are n'^iX) = e^^^''^^^\ Thus, 

(23) ak + 4 - trM^ = e''^'^^^) + e^^'^'^^^) =. 2 cos kW{X). 

Next we prove (|20p. The relation between Cauchy data for the left-to-right scattering solution 
at a; = and x = kl are given by 

au 13k \{ 1 + ^fe \^{tk\(ak{l + rk)+Pk{^{l-rk))=tk 
Ik S, J[i-^{l-r,) ) (^ f^ ; ^ I ^,(l + ,,)+4(^(l-r,)) = ff 

By dividing the second equation by the first one we arrive at 

lk{l + rk) + ^{l-rk) z 
ak{l+rk) + ^{l-rk) ~ bi 
Solving for rfc we obtain 

^ 4-a.-»(&i7fc + ff) 

Using pg|) we get 



sin kW 

sin IB 



Tk 



(«-'^^) + »(^i^mf7 + ^i^/3) 

2cosfc/0 + z(6i^|i^7-^^g^/3) 
This justifies ^. 
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In order to get (PT|) we note that pT]) implies that 



1 + r \ / 1 \ 1 + 7- 1 

= c I ^ 



and therefore, 



f -m+ 
r = ^^^ . 

r- + m+ 



From here and ([T6| it foUows that 



'- - '—^ a+j-P- e''« 



r 



I I e""— a 



bl ' /3 



e''« - (a - ^/3) 



and 



Hence, 



and 



bl e'^^-S b 



^(e*'«-5)-7 e'^'« - (J - zfei7) 



^ + ?T^ 1^ (e''' -S) + l e^'' - (^ + 7^17) 



r 2e''* - (a - -/3) - ((5 + ^6l7) = a + -/3 - (5 - ^617) 
Ol J bl 



(a + ^/?) - (5 - z6i7) (a_^) + ,(A+5i7) 



2e»'9 _ (c, _ ^^) _ (^ + j5^^) 2e»'« - (a + 5) + z(|- - 617) 
(a-J) + i(A+b,7) 
2isin;6' + i(|- -617)' 



where the last equality is a consequence of (TM)) and it implies ((2T|) . 
We prove the third statement of the theorem. From (IMl) it follows 

Using that jrfep + |ifep = 1 we get 

(26) |tfe| :=l-|rfe| = 



Pk \2 ' 



\tk\' 



We use det M^ — akSk — Pklk = 1 and get 



4 + (4 - ak? + (&17. + If )2 ^f^ ((a - J)^ + [b^ + A)2) + , 
From formulas ((25)) and p6| we get 



l^feP 1 ^.. ^2 , .. , A 



(4 - auf + (6i7fe + T^)' 



|t,|2 4V^^ ^^ ^^'^ 61 
and hence, putting k — 1^ 

sin^ /e |ti|2 ^ -^ 
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The analyticity of rfc(A) and r(A) in C+ and their continuity in C+ follow from the explicit 
formulas pO|. (|2T|) . For A S C+lJU^ig„, we have Im6'(A) > and 9{X) is pure imaginary on 
the gaps g„. Furthermore, if Im6'(A) > 0, then 

cos k9(X) 
sm fc6^(A) 

and this justifies the convergence of rfe(A) to r when A G C+ IJ U^]^5„. 

The weak convergence for A G U^]^6„ is a consequence the convergence in the complex half 
plane. 

The proof of Theorem [2] is complete. | 

Note also the following relations: 

|rfe|2 _ sin"^ Ike \ri\^ 



\tk\' = 



aHde 



2 



4cos2fc/0+(6ii^7-^i^/3) 4sin2/0f^S+(5i7-^) 

The last formula follows from \2Q\ by using ([23]). 

Formula p2p implies that the perfect transmission (|ifep = 1) occurs whenever |rip = 
(|tiP = l)orif 

(27) ^^ = 0. 

^ ^ sin^ le 

For e [0, IT /I] equation (P7|) is satisfied when 6*/^ = tott for to = 1, 2, . . . , fc — 1. 

Therefore, in the general case (|rip ^ 0), the transmission probability has fc — 1 peaks with 
|ifcp = 1 in each allowed energy band as 9 increases by ir/l. Since the peaks in the transmission 
probability (or in general in the cross section) are associated with resonances, we expect to find 
k — \ resonances near each allowed energy band. 

On the gaps, the function i9{X) is real valued. Then the transmission probability is given by 

1 



2 



sinh^ lki9 \ri\ 



1 



(28) M\)\ _^_^__^ 

sin\i-'U9 jtiP 

As sinh^;fci6' 7^ for 6* 7^ 0, then in the forbidden zone |ife(A)P ^ 1 unless |iip = 1, |rip = 0. 
Hence there are no resonances below the gaps. 

On the gaps, the reflection coefficient for the half-periodic system r(A) satisfy 

IKA)P = — ^^ ^ 2 = ^ ^ 2 - 1. 

e"^^*^ + e-2»s - 2 + (617 - ^) (a + 5)2 - 4 + (fei7 - §-\ 

In [7], Theorem 5, was shown that 

Lemma 1. // Aq is a degenerate hand edge, i.e. F{Xq) = ±1, F'{Xq) = 0, then the reflection 
coefficient r^ is zero, rfe(Ao) = 0. 

The proof uses the fact that at any degenerate band edge A = Aq the monodromy matrix 
M\g — ±1 and F"(Ao) 7^ 0. This allows to pass to the limit in ((20|) as A ^ Ao. The numerator 
in the right hand side of (l20t vanishes as A — ?► Aq. The denominator converges to ±2/fc, since 
9{Xq) = mr. 



10 
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Thus the medium is transparent for the plane wave with the frequency Aq and we expect to 
find a resonance A G C_ below the degenerated band edge Aq. 

3. The monodromy matrix for Pk 

In this section we give expressions for the elements of the monodromy matrix Mx. 
Let a{x) be 1— periodic function equal to gq for x £ [0, 1) as in ([T|), ([5]): 



ao{x) 



62 ^ for X e [0, 2:2) 
b^'^ for .TG [.X2,l) 



and 



a{x) 



aoix - j) for X G [j,j + 1), < j < fc - 1 
bi^ elsewhere 

for fc > 2, where k is the number of identical cells. The period 1 = 1. The normalized solutions 

ipi{X,x), '(/'2(A,x), 

-dM^)dx^r = A'V'», ^i(A,0) - {b-^d^^l;2){X,0) = 1, ^2(A,0) - {b-^d^^Pi){X,0) = 0, 
satisfy (see [4]) 



V'liA.a;)-^ 2-e:::^-' + 2e 



1 g«A62x _^ 1 e-i-^bs:-; for a; £ [0, X2) 
j^giXb^x _^ ^g-»A6ix for a; e [a;2, 1) 



'(/'2(A,a;) 



ib3.g»A62:r ^ ^g-^A62X for ^ g [0, X2) 



2A' 



2A' 



(^giAbix _,_ j^^-iXb^x for a; g [a;2, 1 

with Aj B, C, D chosen such that "0^ and a{x)dxipj are continuous at X2 

1 



A = — 

4&2 



B = 

C - 
D 



1 

i 

4A 
i 

4A 



(62 + 6i)e*^^^(''^-''i) + (62 - 6i)e"*^^^(''^+''i) 
(62 + bi)e-'^'=^^''^-''''> + (62 - 6i)e'^^^(''^+^i) 
-(62 + &i)e*^"^2(''^"''i) + (62 - 6i)e-'^^^(''^+''i) 
(62 + bi)e-'^''''^''^-'"'> - (62 - 6i)e'^"=^(''^+^i) 



We get the monodromy matrix 

a P 



M.= 



7 S 



V'i(A,l) AV2(A,1) 

i(^a.Vi)(A,l) (^9.V2)(A,1) 



with 



a = -0i(A, 1) = — rr — COS A[5i (1-2:2) +a;252] H — cosA[6i(l - X2) - X2b2\, 



/? = AV'2(A,1) = 



2&2 

&2 +fcl 



262 

sinA[6i(l - a;2) + a:2fo2] -z — sinA[5i(l - X2) - X2b2]^ 



7 = t(i2^2:V'i)(A, 1) = — — — — sinA[6i(l -2:2) + 2;262] „, , sinA[bi(l ~ X2) -2:262], 

A of 201O2 2O1O2 

<5=(^9.V^2)(A,1) ^ 



cosA[fei(l - 2:2) +2:262] Trr — cosA[5i(l - X2) - 2:262] 



26i 



2bi 
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Then 

tT{Mx) = a + S = {p+ l)cosA[6i(l - X2) +2:262] - (p - l)cosA[6i(l - X2) - X2b2], 
where p = -^-g^- The Bloch quasi-momentum 9 = 9{X) satisfy 

2 cos 6* = trAfA- 

The first formula in (|22|) imphes that the one ceh medium is perfectly transparent: |ti(Ao)P = 
1, |ri(Ao)P = 0, if a - (5 = 0, 617 + ^ = 0. We get A = Aq = nm/x2b2. 

Note that there is a resonance A € Res (Pi ) for the one cell operator such that Re A = Aq (see 
equation ((55|) with 63 — bi and xi — 0). 



If An = -^ then we have 
" 2:202 



il-X2)bi 



2F(Ao) = 2cos0(Ao) = ±2 cos ( ; nm ) , if m is even (odd). 

V a;262 / 

Thus in general situation, ^ ^^^' ^ ^ Q, non-rational, Aq is an interior point of a spectral 
band. 

If a;2&2 = (1 ^ 2:2)61 then the Lyapunov function satisfies 

(29) 2F{X) = {p + l) cos(2A62.X2) - (p - 1) 

and we have F(Ao) = 1, F'{Xo) — and F"(Ao) ^ 0. Hence Aq is degenerate band edge. 
The non-degenerate band edge is then given by the equation 

(p + 1) cos(2A62.'z:2) - (p - 1) = -2 ^ A = — ( ± arccos ( | + 27rn 

4. Explicit construction of the resolvent and resonances 

In this section we define the resonances as the poles of the analytic continuation of the resolvent 
R{X) toC_. 

4.1. Representation of the resolvent. In this section we revue some formulas used by Valeria 
Banica in [4 , where she considered the local and global dispersion and the Stricharts inequalities 
for certain one-dimensional Schrodinger and wave equations with step-like coefficients. The 
systems are described by the one-dimensional Schrodinger equation 

{idt + dxa{x)dx)u{t, x) =0 for {t, x) g (0, 00) x 
u{0,x) = uo{x) e L^ 

or by the one-dimensional wave equation 

( {df - dxa{x)dx)v{t, x) =0 for (i, a;) G M x M, 
(31) i v{0,x) =uo{x)eL'^iR), 

[ dtv{0,x) =0 

for a positive step-like function a{x) with a finite number of discontinuities. 
Consider a partition of the real axis 

—00 = xo < 2:1 < 2:2 < . . . < 2:„_i < Xn — 00 

and a step function 

a{x) = b^^ for x € (xi_i,2:j), 
where bi are positive numbers. 

The operator P := -dxa{x)dx defined from {u G H'^{M.), ad^u € H^(R)} to L'^{R) is self- 
adjoint. For ImA > 0, we define R{X)v — {~dxa{x)dx — X^)~^v the resolvent of P. 



i^^) ^ .,(n ^\ _ ...(^\ d T-21 



(32) «(A,..,., . - S C.-«" / #3J^^«, - £ ^V*'-"*, 
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Our choice of the spectral parameter A is related to u in. ^hy X — iuj. We use the expression 
for the resolvent obtained in [4]. 

On each interval (a;i,Xi-|_i) R{X)v is a finite sum of terms 

j{y) e^^^'-y ^ r <y) , 

finite -^f"'' 

where I{xi) is either (—oo^Xi) or (0:^,00) and 

deti:'2(-iA) = {bi + b2)e'^'''^^'-'"\ 
In [4] Banica defines all det £)„ by induction. Let 

detD^HA) == (62 - 6i)e-*^^i(''i+^=). 
We have the following induction relations for n > 3 
(33) 



deti:'„ = e''^''"^"-i 
detA, = 6^*-^''"^"- 



(&„_i - &„)e*^^"-i""-i det Dn-i - {bn-i + &„)e-*^^"-i""-i det A.-i 
(fe„_i - 6„)e-'^^"-^""-^ detZ?„_i - (fe„_i + fe„)e'^'"-^""~^ det A^ 
We define for n > m > 2 






detl^m bm-l + bm 

Then we have for n > 3 

det A^(^A) = (&i + 62)e^^('^^'^)"^n,^2...n-i(&, + 6j+i)e-''^(^^-''^+i)"^(l - dJQ,(-^A)). 

We have induction formula on the Qm's 

(34) g„(-^A) = e^^A^^(^..-^^-0 -^-i+Q-i(-^^) 

1 - «m-l(ym-l(-«A) 



Note that a linear-fractional transform on the unit disc occurs in ([34]) for ImA > 0. 
If ImA < then |e2i^''™(^™-^'"-i)| > l. We use that |d„| < 1 and for any n wc can find e„ > 
such that for every complex A with 

ImA > — e„ 
the estimate 

|Q2(-a)| = |die2*^''^("^-"i)| <1 
holds and gives by induction 

\Qm{-iX)\ < 1, n>m>2. 
Hence (det Z?„(— iA))~^ is uniformly bounded and well defined in this region, which contains the 
real axis. Therefore iXR{X)uo{x) can be analytically continued. The spectral theorem gives 

Lemma 2. The solution of the Schrodinger equation \30^) verifies 



u{t,x)= I e''^'XRiX)uoix) — . 



The solution of the wave equation i31\) verifies 



v{t,x)=l e'*^iXR{X)uo{x) — . 

, ZTT 
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Due to formula ([5^. the resonance spectrum Res (P) consist of zeros of det Z?„(— iA) or equiv- 
alently the zeros of 1 — (iri-iQri-i(— iA) as by ([33]) . equation detZ?„(— iA) = is equivalent 
to 

(3„_i(-iA) = . 

Un-l 

By considerations before Lemma [U for each n> 2 there is e„ > such that all resonances verify 
ImA < — e„. For n — 2 there are no zeros. 

For n = 3 the resonances are solutions of the equation detl?3(— iA) = 0, A = Ai +iX2, Xj G K, 
with constant imaginary part: 



, , 7rm 1 

2b2[X2-Xi) 2b2[xi-X2) 



{h2-h){h2-b^) 



ib2 + b3){bi + b2) 



<0, 



where m = 0, ±2, ±4 . . . , is even, if (&2 — fe3)(&2 — ^i) > or m = ±1, ±3, ±5 . . . , is odd, if 
(&2-53)(&2-6i)<0. 

For 71 > 3 the resonance spectrum can be obtained numerically using the induction relations 

mi- 



4.2. Locally periodic media. Suppose that the profile of a{x) consist of a finite number iden- 
tical elements, obtained by juxtaposing of k unit cells. Outside the interval [0, k) the coefficient 
is constant. We make the following choice: suppose n = 2A: -I- 1 odd, 6i = 63 = 65 = . . . = &„, 
62 = &4 = . . • = bn-ij and put 

si = 0, < a;2 < 1, X3 = 1, d= -^ — -i, A = Ai + 1X2, Ai € R, A2 < 0. 

02 + Oi 

We have X2k+i — k, X2k = k — I + X2- Let 

.. -2 
(36) ao{x) 



62 for X e [0, a;2) 
6j^^ for X e [x2, 1) 



and 



(■37-) ^c ^ ^ f ao(a; - j) for x e [j, j + 1), < j < fc - 1 

^ ' 1 ^r^ elsewhere 

for fc > 2. Here fc is the number of identical cells. 

The function a{x) is called locally periodic or finite periodic on the interval [0, k) with k cells. 
Then, with n = 2fc -I- 1, 

detD„{-tX) - (61 + &2)""'e-^^(''=-''i)'="^n,^2,...,„-i (1 - (-l)Mg,(-7A)) , 
where 

(38) Q2.(-*A) = e2^A..(..-i+..))dct^^ Q2.+i(-*A) = ,2.^,.^^!^ 

det L'2fc det D2k+i 

For fc = 1,2,3,... the resonances are the solutions of the equation 

(39) deti?2fc+i(-iA) = 
or equivalently 

(40) Q2kHX) = 3, 
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1 - aQ2k-2 

Note that if condition ^ is satisfied: 622^2 — bi{\ — X2), then Qm(— *A) as function of Re A is 
periodic with the period TT/b2X2, which imphes the same property for any solutions of equation 
(I40p . The resonance spectrum is then periodic as stated in the Introduction, Theorem [TJ 

5. The poles of analytic continuation of the reflection coefficient rk to C_ 



where 






(41) 


Q2k 


_ ^22Xb2X2 d+Q2k-l 

1 + dQ2k-i 




02 = 


^ ^2l\b2X2^ 



Let rk be reflection coefficient (|24p . In this section we express rk using the iteration formulas 
in section m and extend rfe(A) to A e C_. We show explicitly that the resonances defined in (15^ 
or (Uni) are the poles of r^ , 

— = ^ det2?2fc+i(-iA)^0 ^ 02fc(-iA) = -, 
rk d 

We consider the equation Piu = —dxa{x)dxU — X^u corresponding to the system with one 
unit-cell, a{x) = ao{x), for x G [0, 1) and a ~ l/bf outside [0, 1). Taking solution 

(42) u{x) = Aoe'^''!^ + A^e-*^''^^, a; < and u{x) = ioe*^''^"' + A'^e''^''''' , x > 1, 

the matching conditions imply 

where T is called transmission matrix (see [1]). 

Then coefficients of the solution of the problem with fc— unit cells PkU = ~dxa{x)dxU — X^u 
can be calculated by iteration. Taking solution 

(44) u{x) = yl_ie*^''i^ + A'^^e-'^''''' , a; < and u{x) = Ake"^'"'' + A',^e-'^^''\ x > k, 

the coefficients of the two external regions x < and x > k are related by 

where 

(46) nX) = (^ ' , ^,Jl,,,,^ ) Q'= ^ " 

where 



Q ^il{k+l)\bi J "^ \ f.-il^bi 



Q = ( :^j;[^j ^Jj[^j ) , detg = 1, Qi2(A) = g2i(-A), Qn(-A) = Q22(A) 



is called iteration matrix 

(47) Q{X) = D{X)TiX), D = ^ ^ ^_,,,^ 



^iXbi 



Next we relate the iteration matrix Q and the monodromy matrix M\. 
Relation between the Cauchy data and coefficients in (|44l) is given by 



«(0) \M I \( Ao\_. ( A, 
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and 

"(1) \ M 1 V ^'^'^ \f^0\_rr^fA0 

M,,d(^ j - U -^ J I ^-"^^ J UJ ''-''[ A', 
Using 

u(0) \ / «(1) 

«)(A,0) I ~ 1 (a^xM) 

and definition of T, equation (H51) , we get 

(48) Ma = LDTL-^ = LQL-\ with L ^ [ ]_ _\ 

V bl fcl 



^A I (aa^u)(X,0) I - I (aa^«)(A,l) 



We have 



As 



then 



1.2 _ l2 

gi2(A) = f,^K-^~^^)^^ Jl 1 (e-«A62^2 _ ^i\hr^-i\ ^ 

Q22(A) = e-*-^''i(i-^^) ((r + l)e-'^^^^^ - (r - l)e'^''^^^) . 



detLig = e'^^^^^ ((&2 - 6i)2e*^''^^^ - (62 + foO^e"*^''^^^) 



gJAbi _ g-iAbi 



(49) gi2(A) = —— detZ?3, Q22(A) = -— — dcti^s- 

4O1O2 40i02 

Using the iteration relations between det£'2fc-i, detD2fc-i and detD2fc+i, det£'2fc+i we get 
Lemma 3. We have 

^ Ul?(-A) g<^2^(A) ]' 



where 



kiXbi __ -ikXbi 



L-Hl^L 



Using ([i5|) we have also 

Uo^k + I3klj ~ ^hk + 4^) U^k - Pki-J - ^(7fc - 4^) 

hiak + /3fe^) + %(7fe + 4^) U^k - hi;) + ^(7fe - 4^) 
Using equation p4| we get 

^ (afc-4) + »(bi7fc + t) _ 2Q(y(A) _ q(^,^(A) 

afc + 4 + ^(&i7fc - t ) 2g(^2) (A) Q« (A) ' 

Let first fc = 1. Using equation p9|) we get 



—iA(l — 2:2)61 pi A6ia;2 

g2i(A) = -e-*2A(i-.2)biQ^2 = — jt^t:^ det^3(A) 
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Hence 

g-2iA6i(l-X2)giAbi(l-X2)giAbia;2 jjg^ f) 
^^ ~ g-jA6i det Dg ~ 

detDa 
For fc > 1 we use 

and get 

Qf) = _g-.2A(l-.2)b.QW ^ _g-.2A(l-.2)b.^_^det^2.+l. 

(46162) 

Hence, by (jM)) . 

'^ qW e-*^-^''idetZ?2fe+i detZ?2fc+i 

We have 

(50) rfc = -e 1^ "^Q2/c+i -*A = .„ , ... . 

Thus we have proved that the poles of r^ on C_ coincide with the sohitions of (|40|) : 

1 



l-dQ2fe(-JA) = 4^ Q2fcHA) = ^. 



Hence, we have 



Proposition 1. The reflection coefficient rj.(A) continuous to C_ with the poles at the resonance 
spectrum Res (Pk)- 

6. Convergence of the resonances to the real axis as fc ^ 00 

In this section we prove that the resonances spectrum Res(Pfc) converges to the real axes as 
fc — >■ 00. 

First we note that the function 9{X) — arccos ((a(A) + <^(A))/2) has analytic continuation 
onto the domain C \ IJ^i 5" ^y ^^^ formula 9{\) = 9{\) and Im0(A) < for ImA < 0. For 
A € U^i 9n ■we set 9{X) = 9{X — iO) and 9 is pure imaginary there. 

By Proposition [1] the function r/c(A) is analytic in A G C_ \ Res (Pfc), and the formula 

^ (a-J) + »(6i7+£) 

'''k[X) ^ , a, -.-, cos kW(\) , ■/, 8\ 

2sm?g(A) ^.„^^g(\; +z(6i7^^) 

extends to C- \ Res {Pk), where we have Im0(A) < 0. 
Furthermore, if Im0(A) < 0, then 

cosfc0(A) 

. ,n,x\ ^ « as fc ^ CX3, 

smfcy(A) 
and we have 



2sin/0(A) + (6i7-^) 



rk ^ /. ,'. \ "v. = f{X). 



Note that for A e R, r{X) = r(-A). 
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The limit extends also to A e Uj^j^^n, where 9{X) is pure imaginary. 

Let gfe(A), g{X) denote the denominators of r^ respectively r for A G C_. Then 

5fc(A) ^ g{X)^2smW{X)+(hj-^^ 

uniformly on any compact subsets of C_ and the limiting function g is analytic on C_. 

For any k all zeros of ^^ have negative imaginary part: Res (Pfc) E C_. Then, by the Hurwitz's 
theorem, the zeros of gk can only converge to the real axis. | 

7. The limit limfc^oo rk as a fixed point of a sequence of linear-fractional 

AUTOMORPHISMS OF THE UNIT DISK 



Equation (|50p together with iteration relations ([50| shows that for A G M real and for all 
k = 1, 2, 3, . . . the reflection coefhcients for the fc+1 and k cells media are related by r^+i — fx{rk), 
where f\ is linear-fraction automorphism of the unit disk. Hence we get 

(51) rfc = /f(ri):= /a o /;,... /;,(ri). 

k iterates 
Here ri is the reflection coefficient for Pi, 

^^ ~ ~ 1 - (pg2i\b2X2 

and for A G M we have 



|2 



A(f 



In this section we will consider the limit of the sequence /| as k ^^ oo. 

First we recall some well-known facts on the convergence behavior of a sequence {/''^'} when 
/ is general linear- fractional automorphism / of the unit disc I? = {zgC; |z|<1}, 

/(z) = -^, heD, zgC\{1/6}. 

We refer to the paper of Burckel [5]) for the details. 

In general situation (6 7^ 0), / has two fixed points zi, Z2- There are three cases to consider. 

Hyperbolic: / has two (distinct) fixed points on dD 

\zi\ = 122! = 1, Zl 7^ Z2. 

In this case the sequence {/''^'} converges uniformly on compact subsets in D to one of 
these points. 
Parabolic: / has one (double) fixed point on dD 

zi = Z2, with \zi\ — 1. 

In this case the sequence {/f^'} converges uniformly on compact subsets in D to this 
fixed point. 
Elliptic: / has two fixed points: one fixed point zi & D and one fixed point Z2 — 1/z ^ D. 
In this case either / is periodic in the sense that /'"l = / for some n, or the orbit 
{/' '; n G N} is dense in the compact group of all conformal automorphisms of D which 

fix Zi. 
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We apply these results to f = f\. In order to simplify the formulas we suppose 

(52) b2X2^h{l~X2). 

Then we have the following expression for f\ : 



-d + r] 



d-\-rjZ 



h{z) = ^iT^, where r; := e^^^^^^^ 

In general situation ?] 7^ 1, equation f\{z) = z have solutions 



. . _ -(l + ??)± yj {I + -qf - AdPrj _ - cos(Ab2a:2) ± Vcos2(Ab2a:2) - (P 

Note that 

(54) if d^ < cos^(A62a;2) < 1 then \zi\ < 1, |z2| > 1 and zi-z^^l 

(55) if cos2(A62a;2) < d^ then |zi| = |z2| = 1, 

and zi = Z2 if cos^(A62a:^2) = o^- 

Note that if ((52l) is satisfied then the Lyapunov function is given by (|29l) . -F(A) = (p + 
l)cos^(A62a;2) — p. The inner points of the allowed bands (JJ^^ bn satisfy |F(A)p < 1 which is 
equivalent to rf^ < cos^(A62a^2) < 1- On the spectral gaps IJ^i 9n we have cos^(A62a;2) < d"^. The 
non-generated edge point of a band is given by cos^(A62X2) — d"^. The degenerated edge point Aq 
satisfy cos^(A&2a;2) = 1 ^ Aq = Trm/b2X2- 

Hence we get the following version of the forth result in Theorem [2| 

Proposition 2. If X (z Un=i 5" then the linear fractional automorphism f\ has two fixed points 
zi, Z2 of hyperbolic type: \zi\ = \z2\ = 1. The sequence rk converges to r as k -^ 00, where 

is either z\ or Z2. 

//A is a non- degenerated band edge point: F{X) — ±1 and F'{X) ^ 0, then f\ has one (double) 
fixed point z\ = Z2 and \z\\ = 1. 

We have 

r = zi = lim rk = lim /| (ri). 

// A = Aq is degenerate band edge point F{Xo) = ±1, F'{Xq) = (degenerate band edge) then 
r-fc(A) ^0 for all k^ 1,2,3.... 



The Proposition is still valid if condition ([52]) is not imposed. 

Appendix A. Numerical calculations 

In this section we present some examples of the resonance spectrum. The resonances are 
solutions of the equations _D2fc+i(— «A) = 0, where _D„, n = 2fc + l, are defined iteratively by p3l) . 
Section |4l The zeros of I?2fc+i(~*A) are calculated numerically by using the Newton procedure. 
Using Matlab we plot the resonance spectrum for the number of identical cells k — 3,4, 5. In the 
same figure we show the band spectrum for the corresponding periodic problem satisfying (|4|. 
The small circles on the real axis marks the position of Aq = nm / {x2b2) when the one-cell system 
is perfectly transparent: |ii(Ao)P = 1, |ri(AoP = 0. 

On Figure ^ condition (jB)) is satisfied: 61 = 1, 62 = 4, X2 — 0.2. 

On Figure ([2]) condition (jB]) is not satisfied: &i = 1, &2 = 3.8, X2 = 0.2. 

On Figure ([3]) condition ^ is not satisfied: bi = 3.8, 62 ~ 1, X2 = 0.8. 
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